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Abstract 

The correlation properties of the cold system consisting of the electron-hole plasma interacting with the exciton gas are 
analyzed. It is shown that the homogeneous state of the system is unstable and in the stationary state the densities of the 
electron-hole plasma and exciton gas are modulated. 
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The problem of the existence of the Bose condensed state of the exciton gas in semiconductors is the subject of the 
experimental and theoretical investigations for a long time. However, there is no undoubted proof of the existence of 
this state until now. 

Recently, the properties of the cold exciton gas are investigated intensively, both experimentally and theoretically, 
in double 2D quantum wells, besides, electrons are localized on one plane of these wells and the holes are localized 
on the other plane. The state of the cold exciton system being observed in the experiments of different experimental 
groups [H, 0], 0, 0, [1], [13], ill, [1], [1], is nonhomogeneous, and the density of the system in these experiments 
is modulated in space. 

The theory of this phenomenon based on the Turing kinetic mechanism of the instability [l^ have been proposed in 
works 0, 0- In this theory the system of nonlinear diffusion equations with sources of particles for electrons, holes 
and excitons are proposed and solved. The solution of these equations demonstrates the state with the periodically 
modulated density of the exciton system only if the constant describing the decay of excitons depends on their density. 
However, in the low density approximation this dependence should be neglected. Moreover, this mechanism has the 
classical character and does not describe the systems with the quantum coherency, but for the sufficiently small 
temperatures one can suppose the existence of the coherency in the exciton system , , ■ 

Another theory, based on the supposition of the existence of the attraction between excitons at large distances and 
the formation of the liquid phase in the low-density exciton system due to this attraction, has been proposed in the 
work The supposition of the existence of the liquid phase in the low-density exciton system is baseless. The 

existence of the attractive part of the interparticle interaction potential in 2D Bose systems results in the collapse of 
the system to the densities of the order of the squared inverse radius of the repulsive core [? ]. However, for these 
densities excitons can not exist as individual particles, they are destroyed and form the electron-hole liquid. Note that 
in the usual 3D case, as it was shown in the work [IS], the scattering amplitude for the isotropic excitons is positive 
and is equal to the value a = -S^ob, where as is the effective Bohr radius, and just the scattering amplitude, but not 
the bare potential, determines the properties of the system in the case of small densities. The effective Bohr radius 
has the form an = — r — vi\ where h is the Planck constant; m* = is the reduced mass, mg u are electron (e) 

and hole (h) masses, which are supposed to be isotropic and of the same order ^ mh', e* is the effective charge 

/ 2\l/2 

in the semiconductor ^* = [j^j ! where e is the bare electron charge and xq is the static dielectric constant of the 
semiconductor. As the result, the 3D exciton system of small density can exist only in the gas phase. The repulsion for 
the interaction between excitons in 2D system, where electrons and holes are localized on different planes, is sharply 
defined than in the usual 3D case due to the geometrical constraint. In this case the formation of the exciton liquid 
phase is possible only for the densities n^x > ^/o-ex^ where a^x is the size of the exciton. The size of the exciton is 
equal to the effective Bohr radius aex = in the case of the small distance I between the planes of the quantum wells 
I << Ob, and aex » as for the large distance between planes when I » as- If the exciton density is sufficiently 
large n^x > ^/o,ex the excitons can not be considered as the individual particles. In this case the electron-hole system 
represents the strongly interacting electron-hole liquid. 

In the present work the Coulomb correlations in the cold electron-hole plasma as well as the quantum effects and 
the coherency of the exciton subsystem in quantum wells are taken into account. We suppose that the external source 
is stationary, has the frequency larger than the semiconductor gap Ag and creates electrons and holes being localized 
on different planes. Moreover, the region of the action of the electromagnetic field is supposed to be sufficiently small. 
In this case the electron-hole plasma being created by the external source propagates along 2D planes, besides, the 
electrons propagate along one plane and the holes along the other plane. During this propagation the electron-hole 
plasma loses the energy giving away the energy to the phonon system and comes to the equilibrium state with the 
small temperature T. Note that the parameters of the stationary state, such that the densities of the electron-hole 
plasma n and the exciton g defined by the equilibrium between incoming rate of the electrons, holes and 

excitons and outgoing rate of these particles. 

Thus, the model being analyzed represents the quasi-equilibrium electron-hole plasma from which the exciton gas is 
created as the result of the formation of the electron-hole bound states. In the present work we consider the properties 
of the stationary state of this system, namely, the properties of the electron- hole plasma of the small density na^g << 1 
interacting with the exciton gas of small density n^xag << 1. 

The small density of the electron-hole plasma and the exciton gas are necessary if we suppose the existence both 
electron- hole plasma and the exciton gas. In the opposite case of the large density na^ >> 1 the excitons are 
destroyed and can not exist as the individual particles. This case is not under consideration in this work. Note that 
the experiments Q demonstrate the existence of the excitons, and this means the existence of small density 

of the electron-hole plasma in these experiments. 

We analyze the correlation properties of the electron-hole plasma of the small density and the interaction of this 
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plasma with the exciton gas. The energy of the ground state of the electron-hole plasma has a minimum value as the 
function of the density at some density uq such that uqu^ ~ 1. As a result, the electron-hole plasma tends to form 
the liquid electron-hole drops and the homogeneous state of the system becomes unstable. The inhomogeneous state 
of the electron-hole plasma creates the periodic mean field which acts on the exciton gas and makes the density of 
the exciton system modulated in space. 

We suppose that the system consisting of the electron-hole plasma interacting with the exciton gas is in the 
stationary state with some small temperature T obeying the inequality T << e^''*' << Ejf^^; where e^''*'' are the 
Fermi energies of the electrons and the holes correspondingly Sp''^"' = ] the Fermi momentum in 2D case is 

1/2 

expressed by the equality pf = (27rn) ; the average density n of the electron-hole plasma is equal to the density of 
each component n = Ug = rih] e\^^^ is the exciton bound energy 

= Below we use the system of units in 

which h = m* = e* = 1. 



EFFECTIVE ACTION 



To analyze the properties of the exciton system interacting with the electron-hole plasma the effective action for the 
slow fields is derived. In this paper the properties of the stationary state of this system is under consideration, but in 
future the kinetics of the formation of this state is supposed to be analyzed. For this reason the Keldysh-Schwinger 
technique for the nonequilibrium processes 17 1, [l9| in the functional integral formulation (see (20| . for example) is 
used for the derivation of the effective action of the slow fields. The generation functional for the considering system 
can be written as 



Z = J Y^DtpaDtl;^ exp (^iSe-h + i<^dt(fr {Jai^a + "(paJa, 



(1) 



The action of the electron-hole plasma in the electron-hole representation has the form 



^Tpaix) idt+ tia- £a{y^ Ip a {x) 5 [x - x') 



S^h = (fdrxdrx < 



Q.—e,h 

. 1 ^ (-1)"+^ {x) ^« (x) Uc^p (x') ^Pp ix') S (t - t') 

a, 13 



(2) 



here x = (t, T^) denotes the collection of time t and space coordinates in 2D; "ipe and iph are fermion (Grassmann) 
fields of electrons and holes; £e ("P*) ^^'^ ^'^'^ dispersive laws of the electron and hole bands £a ( P*) = 
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2^1— ~p* , where Ag is the semiconductor gap; ~j} — —iV is the momentum operator; Uc — Uc C?^ — ~r^') = 
is the Coulomb interaction. In the Keldysh-Schwinger technique the time variable t changes along 



Eo\r - r '\ \r - r'\ 

the double time contour with return. 

The effective action for the electron-hole plasma with the formation of the exciton system have been obtained in 
the functional- integral technique in (23j , 24 1 . In this paper we give more simple derivation of the effective action and 
introduce some simplification transforming to the density functional for the Fermi subsystem. Integrating over the 
rapid ipa-&e\ds, changing on the scale much smaller than the average distance between the charge carriers in 
the ladder approximation we can obtain the action for the smooth fields in the form 



■0Q (idt - ^a) IpaS (x - x') - 



d^xd^x' 



a—e,h 



\ (V'eV'e - V'ftV'/i) Uc (V'eV'e " i^h^h) S {t - t') - 



aJ3—e,h 



(3) 



here the fields ipa are the smooth Fermi fields. The vertex Tap is the sum of the ladder diagrams. The internal 
Green functions of these diagrams have the large momentums, much larger than pp, and the external lines of these 
diagrams have the small momentums, much smaller than the internal lines. The value of the boundary momentum A 
separating the rapid and the smooth fields obeys the inequality ^ >> A » pp. The vertex F can be considered as 
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independent of the external momentums or frequencies due to the large values of the niomentums corresponding to 
the internal lines of this vertex, on the assumption that this vertex does not contain the pole. The vertexes T^e and 
Fhh do not contain the pole parts, but Teh has the pole part corresponding to the formation of the bound state of 
the electron and the hole, i.e. the exciton. 



2m* 



i/"" ( ) C ("g ) r (f^ V'y ( ) V'y Cq) 



E- 



Ei2'^+id 



E 



E^ + id 



Uc (y - t ') I (4) 



The wave functions (^k^ 
electron- hole pair, -ip^ ^ fc ^ 



are the relative motion wave functions of the discreet part of the spectrum of the 
are the wave functions of the continuous part of the spectrum of the electron-hole pair; 



the energies Ee'x are the spectrum of the bound energies of the exciton; here we denote E = Q, + He-h — P"^ /2M, 
where He-h = Pp/2m*. The incoming "p e, and outgoing "p g, "p^ momentums of the ladder diagrams can be 
represented in the form 



'p: 



M 



Ph- ^P k 



Ph 



M 



k' 



where the momentum P = 'p^ + ~p ^ is the exciton momentum, i.e. the total momentum of the electron-hole pair 
forming the bound state. The wave functions ipn ("^) and (T*) of the discreet and continuous parts of the internal 
exciton spectrum obey the Schredinger equation 



2m* I f 

V2 1 



^ j (^) = E^Hn (^) 



2m* 



^1? i'r) = Eex Cp) tp^ {'r) 



The normalized solution of this equation in 2D for the lower energy state can be written in the momentum repre- 
sentation as 



V'o {~p) 



(1 + ^2) 
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The normalization of this wave function has the form J (f^^'^^o ("P ) = 1- 

We suppose that the energy E is near the lower bound energy of the exciton e'^^ . Separating the pole and the 
non-pole terms we can write the vertex Teh as 

(c) . 

where the term F^^ is the non-pole part of the vertex Teh, which can be considered as a constant of the order of 

unity T^^f^ ~ 1; the collection of the summarized frequency and the summarized momentum is P = (n, . Due 

to the smallness of the incoming and outgoing momentums | 'pe,h \, \ ~p'eh l<< = 1 ^^e vertex T^h can be 

represented in the form 



O + lle-h - ^ - Eil' + id 
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where the constant F in the case of the small incoming and outgoing momentums can be represented as 



= (e(°)Vo(0))' 



(c) 

Near the pole of the vertex Teh the term F^^' ~ 1 can be neglected. The vertex Teh can be decoupled by the 
introduction of the virtual exciton field b, and the action of the system can be written in the form 



where 



^e. = ^-^le-h + E<f) (8) 

Eb is the bound energy of the internal exciton ground state. As it will be seen lower, the chemical potential for 
excitons is connected with He-h as 

He.=He-h-E<^J-<T^'^l> (9) 

at that, the value < T^*^ > can be written as 

<TW>=GW(a.-a;')G(°)(x-a;')= / ^} (10) 

\-p\>{^ 

It is convenient to transfer from the Fermi fields of electron-hole plasma to the density variables, i.e. the density 
functional with respect to the Fermi fields. This transition can be accomplished by the introduction of the functional 
(5-function to the statistical sum 



Z = j Di,c.D^„DbDbexp{iSe,h,ex) = 

= j D^^D^^DbDlDuo, exp {iSe,h,ex) J{5 [V^„^ 



Using the Fourier representation for the functional J-function 



Jj5 [^aV'a - ria] = j DVa exp d^xVa {x) [V'aV'a " ' 



we can rewrite the action in the form 

-2_^Vocna - ^{rie-nh) Uc (rie - n/j) + b {idt - ^ex) b 
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where 



Z = y D^aD^J^DbDbDriaDVa exp {iSe.h.ex [i'a,b, Va,na]) 

It is convenient to introduce new variables instead of Ua and Va, namely, 



(12) 



n + Sn 



Ve^V + SV: 



n ~ Sn 



(13) 



Vh^V~dV 



Note that the fields Sn and SV correspond to the fluctuating violation of the electro-neutrality in the electron-hole 
plasma. The integrals over the fields tpa and 6n in Eq. (fT2|l are Gauss integrals and can be calculated. As a result we 
obtain the action in the form 



Sex [b,Va,na\ = -iSpln 



idt ~ -_iv + SV) b 

b idt -^h-{V- dV) 

+ fd^x^Vn + ^5VU(^^SV + b{idt- ^ex)b 



(14) 



The supposition of the low exciton density gives the possibility to expand in powers of the field b. Integrating over 
the fields SV and neglecting the powers higher than 6** we obtain 



Sex [b,Va,na] = ~iSpln 



- 



idt-^e-V 

idt-Ch~V 

2 ~> ( 



(15) 



b 
b 



G(") [V] 



^Sp(- 



b 
b 



G(o) [V] 



Ecorr [V] + <j)d^x {Vri + b [idt ~ ^ex] b} 



The functional Ecorr [V] is the correlation energy. This functional is represented as the sum of all simply connected 
closed diagrams dressed by the internal Coulomb interaction lines and the lines of the external field V. In the case of 
constant V the correlation energy Ecorr [V] can be represented in the well-known form 



Ecorr [^] 



-i-Sp In 
2 



1 - n [y] Uc 



A f dwifk Ucik)Il[V] 



{2Tr)^ 1 - gUc (k) n [V] 



(16) 



where 11 [V] is the total polarization operator. The simple transformations of the second and third terms in ([T 
give the action Sex in the form 



Sex [V,l 



-iSpln 



idt-^e-V 

idt~^h-V 

1 



Ec, 



(17) 



+ (pd-^x^Vn + b [idt - Cex] b - -'jex {bb) - 7ex-e,/i {bb) n 



The value < r|*J] > has the form (jlOp . The coupling constant "/ex is the exciton-exciton scattering amplitude; 
and ^ex-e,h = lex-e + lex-h whcrc ^ex-ei lex-h ^rc the excitou-electron and exciton-hole scattering amplitudes, 
respectively. 

The first order expansion of the first term in Eq. (jl7[) over the field V gives the change Vn — > V (n— < n >), 
where < n > is the average density of the electron-hole plasma in stationary state <n>= — where Vol is 
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the volume of the system. The second order expansion of the first term in Eq. pT|) in series of the field V and the 
integration over this field gives 

SeAl>,n]=L^J _[Hn-<n»{TiiorUn-<n>)-E...An]\+_ 1 ^^^^ 

J [ +b [idt + Hex - f^e-h - £,ex] b - \^iex (bb) - Jex-eji (bb) 71 J 

where H^"^ is the polarization operator of the form If'^) {x — y) = —I'^^Ga ^ {x — y) gL"' {y ~ x). Expending in 

series of the field V we suppose the smallness of this field. This supposition is correct when the fluctuations of the 
density n are small compared with < n >. This supposition is correct for the beginning of the density fluctuations 
growth. The derivation of the effective action in the supposition of the smooth variation of the density n at the scales 
lager than the average distance between particles without supposition of the smallness of the density fluctuations will 
be considered in the next paper. 



INSTABILITY OF THE SYSTEM AND THE DENSITY MODULATION OF THE ELECTRON-HOLE 
PLASMA 



Using the effective action for the interacting electron- hole plasma and the exciton gas Sex [b, ?^] Eq. (fT8|) we show that 
the homogeneous state of this system is unstable, and due to the Coulomb correlations the low density electron-hole 
plasma should form the inhomogeneous state. The formation of the inhomogeneous state of the low density electron- 
hole plasma is the result of the tendency of the low density plasma to collapse to the liquid electron-hole drops. The 
interaction of this plasma with the exciton gas results in the existence of the self-consistent nonhomogeneous periodic 
field acting on the exciton system as the periodic external field. This periodic field results in the density modulation 
of the exciton system. 

To analyze the effective action for simplicity we put m = rUe — riih — 1- The polarization operator II*-'^-' ^w, k ^ for 
I k \«pp IS calculated as 



27T LO' 



The correlation energy Ecorr [n] for small densities na^ «1 can be calculated as 



for Lu >> kVp 



(19) 



Er 



(20) 



where the constant A is estimated as A ^ 1. The main contribution to the integral for Eeorr \p\ (HSl) is given by 
the momentums /cq ^ n^^^ >> pf ^ n^^^ and the frequencies luq ^ k^. Note that the Hartree-Fock energy for this 
system is estimated as 



Eh-f [n]^ - - 



PF 



(21) 



Thus, the contribution of the Hartree-Fock energy can be neglected compared with the correlation energy 
Eh-f (n) << Ecorr ("-) for na^ << 1. 

The substitution of the polarization operator n'°) (^uj, k ^ Eg. p^ for uj >> kVp to Eq. pS)) gives 



Sex [n] 



(An) 



(An) - Ecorr (n) 



[ +b [idt - £,ex] b - \^cx {bb) - jex-eji {bb) ■ 



where An = n— < n > and expanding in An the functional Ecorr [n] we obtain 
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Sex [^n] — (pd^x 



{An) 



2-K (kV„Y 

1 /-x2 

2 Tea: 



ihh) - ^ex-e,h {bb) <n> 
In these expressions the frequency uj should be considered as lo — 



-lex-eJi {bb) An 



(22) 



E' 

corr 

E" 

corr 



idt- The term hnear over An, namely, 
(< n >) An, can be omitted if the total number of particles is fixed, i.e., / An = 0. The substitution of 
(< n >) (HOI) to dig) gives 



Sex [Sn] ^ fd 



At 



' 2^ex 



|A < n >-2/3 



An + b [idf - ^ex] b+ 
lex-eM (bb) <n> -jex-eji (bb) An 



(23) 



From the first term in Eq. it can be seen that the system is unstable relative to the phonon oscillations 

of the electron-hole plasma density. The sound velocity of the density oscillations has the imaginary value | c p = 
< n >~2/3 Vp) . Note that the sound character of the considering density fluctuations is not surprising, since 
these fluctuations do not violate the electro-neutrality of the electron-hole plasma. This instability results in the 
growth of the density fluctuations of the plasma with the increment of the growth 



/A \ 

(:=(^A\ <n>-^^^ kVF k <n>^/^» kVp 

This increment of the growth corresponds to the frequency which is much larger than kVp, C >> kVp. Due to 
this inequality we can use the asymptotic of the polarization operator XI^"^ {uj, k) ^ '-^^a being correct for the large 
frequencies a; >> kVp- 

As the result, the stationary state of the electron-hole plasma should be nonhomogeneous, besides, the modulation 
of the density has the periodic character. 

The modulated density of the electron-hole plasma results in the existence of the self-consistent field acting on the 
exciton system. This field has the form ((23|l 



V{V) =^ex-e,hAn{~r) 

Note that this external for the exciton system potential plays the essential role, if the amplitude of this potential is 
larger than the chemical potential of the exciton gas. This condition can be obeyed for the sufficiently large density 
of the electron- hole plasma, i.e., for the sufficiently large pump. 

At the conclusion we estimate the radius of the character non-homogeneity of the electron-hole plasma from the 
kinetic consideration [i^] . We suppose that the kinetics of the non-homogeneity formation is analogous to the kinetics 
of the electron-hole liquid drop formation or the phase immiscibility and can be described by the equation 



dt (nR^no) = 2TTRVFn + 2TTRVexnex - ttR^— (24) 

The first term in the right hand side of the equation (|24p is the incoming term and the second one is the outgoing 
term; Vpn is the current of Fermi particles incoming into the drop of the radius R; Vexnex is the current of excitons 
incoming into the drop, Vex is the average velocity of excitons which can be written as Vex ~ Vl^^iE; ^ jg iii^ density 
of the electron-hole plasma in the gas phase. Hex is the density of excitons, and ng is the density of the electron-hole 
liquid within the drop; r is the lifetime of the carriers. There are no vapor terms in Ea. (|24p since T=0, and particles 
can not overcome work function from the liquid drop. The equilibrium value of the radius R can be found from 
Eq.(l24l) as 



i?o = 2^^" + ^-"-r (25) 
no 

Note that the radiuses larger than Rq are not profitable due to the different dependence on R of the incoming and 
outgoing terms - the first one is the surface term and the second one has the volume character. From the point of 
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view of the phase immiscibihty the hquid phase is the electron-hole liquid drops and the gas phase is the mixture of 
the exciton gas and the electron-hole plasma of the small density. 

Finally, the set of drops forms the periodic structure. We suppose that this is the consequence of the homogeneity 
of the system in average. This periodic structure forms the external periodic field for the exciton gas. If this Bose 
gas could be considered as the ideal, it should be localized at the minima of the periodic external potential. 
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